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1.  Introduction 

There  is  significant  research  interest  in  developing  new  materi¬ 
als  and  structures  with  improved  survivability  when  subjected  to 
high  strain  rate  loadings  associated  with  ballistic  impact  and/or 
blast.  These  new  material  systems  may  be  created  with  sub-struc¬ 
tures  and  components  of  different  length  scales,  the  combined 
effect  of  which  is  to  better  dissipate  the  energy  of  the  blast  or 
impact.  Such  multi-scale  systems  require  more  advanced  models 
of  the  kind  discussed  in  this  work,  which  are  often  referred  to  as 
micropolar  or  Cosserat  models,  e.g.  Mindlin  (1964)  and  Eringen 
(1999).  Models  of  this  type  possess  additional  kinematic  degrees 
of  freedom,  e.g.,  a  point  in  such  a  continuum  may  be  capable  of 
both  translation  and  rotation.  With  additional  degrees  of  freedom 
come  additional  constitutive  properties  and  stress  terms. 

Improving  our  understanding  of  failure  mechanisms  and  mak¬ 
ing  use  of  this  knowledge  in  developing  improved  predictive  mod¬ 
els  is  another  aspect  of  ongoing  research.  For  example,  regions  of 
armor-grade  ceramics  are  often  pulverized  as  a  result  of  impact. 
The  modeling  of  these  post-failure  comminution  zones  may  be 
improved  by  viewing  such  zones  as  regions  of  granular  media, 
another  field  of  study  where  micropolar  models  allowing  rotations 
of  individual  particles  are  often  employed.  Just  as  a  cohesive  zone 
has  been  used  to  model  damage  in  an  elastic  strip  by  Gazonas  and 
Allen  (2003),  so  also  a  micropolar  layer  could  be  used  to  describe 
the  failed  region  in  a  ceramic  subjected  to  secondary  impacts. 

Various  experiments  have  been  developed  to  characterize 
material  properties  by  subjecting  test  specimens  to  impact  loads 
and  measuring  resulting  deformations  and  propagating  stress 
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waves,  such  as  the  split  Hopkinson  pressure  bar  and  the  plate 
impact  test.  The  plate  impact  test,  for  example,  effects  a  state  of 
uniaxial  strain  well  suited  to  a  one-dimensional  analysis.  The 
analytical  solutions  to  simplified,  one-dimensional  models  may 
give  insight  into  the  behavior  of  materials  subjected  to  such  tests, 
suggest  desirable  constitutive  properties  or  design  parameters,  and 
allow  for  the  determination  of  constitutive  properties  from  exper¬ 
imental  results. 

Although  this  work  uses  a  one-dimensional,  linear,  anisotropic 
micropolar  model  to  analyze  the  effect  of  instantaneously  applied 
loads  or  impacts  leading  to  transient  wave  propagation,  there  are 
many  other  works  available  in  the  literature  that  consider  related 
problems.  The  governing  equations  developed  in  Section  2  are  sim¬ 
ilar  to  those  encountered  in  the  study  of  helical  springs  by  Jiang 
et  al.  (1991)  and  twisted  ropes  or  cables  by  Samras  et  al.  (1974), 
Ostoja-Starzewski  (2002),  and  Shahsavari  and  Ostoja-Starzewski 
(2005),  namely  a  system  of  two  coupled  partial  differential  equa¬ 
tions  (PDEs).  Raoof  et  al.  (1994)  modeled  an  impact  by  applying 
a  step  load  to  a  spiral  strand.  In  much  of  the  literature,  harmonic 
solutions  are  considered  in  the  study  of  wave  propagation  in 
micropolar  media.  For  example,  harmonic  solutions  were  obtained 
for  wave  propagation  along  a  composite  wire  rope  using  coupled 
PDEs  by  Martin  and  Berger  (2002)  similar  to  those  used  in  the  pres¬ 
ent  work.  Krishnaswamy  and  Batra  (1998)  examined  wave  propa¬ 
gation  in  a  linear,  infinite  Cosserat  rod  with  two  directors  and 
examined  the  effects  of  dispersion  while  considering  harmonic 
motion.  Finally,  although  the  majority  of  the  literature  considers 
the  linear  problem,  there  has  been  some  recent  work  studying 
wave  propagation  in  non-linear  systems  by  Porubov  and  Pastrone 
(2004)  and  Pastrone  (2005). 

A  more  general  consequence  of  the  coupled  PDEs  used  in  the 
present  work  is  the  existence  of  multiple  waves.  Applying  a  stress 
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pulse  at  one  end  of  a  one-dimensional  micropolar  bar  will  lead  to 
two  axial  stress  waves  with  different  wave  speeds  traveling  down 
the  bar.  Some  other  examples  from  physics  of  this  phenomenon  in¬ 
clude  birefringence  (the  decomposition  of  a  single  ray  of  light  into 
two  waves)  and  wave  propagation  in  poroelastic  columns,  see 
Schanz  and  Cheng  (2000,  2001).  This  latter  example  results  in  a 
system  of  two  coupled  PDEs  for  the  one-dimensional  case  in  terms 
of  displacement  and  pore  pressure  based  on  Biot’s  theory  of  poro- 
elasticity.  Two  compressional  waves  were  then  observed  in  the 
resulting  one-dimensional  model  of  a  poroelastic  column.  Plona 
(1980)  also  makes  reference  to  Biot’s  theory  when  he  discusses 
his  experimental  observations  of  a  second  compressional  wave  in 
a  porous  medium.  It  is  hoped  that  the  work  presented  in  this  paper 
may  benefit  applications  with  similarly  coupled  systems  of  govern¬ 
ing  differential  equations. 

In  Section  2  of  this  paper,  the  mixed  initial-boundary  value 
problem  for  two  different  micropolar  systems  is  formulated.  The 
first  system  is  a  two-layer  bar  rigidly  fixed  at  one  end  with  a  stress 
applied  at  the  other  end.  The  second  system  consists  of  two  differ¬ 
ent  micropolar  bodies,  a  flyer  and  a  target,  impacting  one  another. 
An  equivalent  discrete  model  of  a  rigidly  fixed  bar  is  included  to 
provide  additional  physical  insight  by  relating  the  model  parame¬ 
ters  from  the  continuum  model  to  the  discrete  model.  The  solution 
methods  used  for  these  problems  are  described  in  Section  3  and  in¬ 
clude  both  a  Laplace  transform  approach  and  a  D’Alembert  ap¬ 
proach  to  solving  the  boundary  value  problems.  Section  4 
contains  five  different  examples  including:  an  analysis  of  a  discrete 
system,  a  comparison  between  the  discrete  and  the  continuous 
systems,  a  number  of  impact  examples,  and  a  study  on  reflection 
and  transmission  coefficients  due  to  impedance  mismatch  be¬ 
tween  different  materials.  Finally,  Section  5  summarizes  the  results 
of  this  work  and  suggests  some  ways  these  results  may  be  used 
and  expanded  upon  in  future  studies. 

2.  Formulation  of  the  governing  equations 

The  general  form  of  the  mixed  boundary-initial  value  problem 
describing  the  behavior  of  a  linear,  one-dimensional,  anisotropic, 
micropolar  body  is  presented  in  Section  2.1.  Two  particular  config¬ 
urations,  a  two-layer  bar  and  an  impact  problem,  are  formulated  in 
Sections  2.2  and  2.3.  Finally,  details  of  an  equivalent  discrete  model 
are  presented  in  Section  2.4.  Although  the  main  emphasis  is  the 
analysis  of  continuous  problems,  the  discrete  problem  is  included 
to  provide  additional  physical  insight  into  the  nature  of  micropolar 
media. 


For  the  one-dimensional  model  under  consideration,  the  posi¬ 
tion  along  the  body's  length  is  parameterized  by  x.  The  strain  mea¬ 
sures  for  the  normal  strain,  e,  and  the  micropolar  strain,  y,  are 


(2.2) 


As  is  often  done  in  micropolar  models,  <p  is  considered  to  be  a  rota¬ 
tion  and  u  is  an  axial  deformation1  as  shown  in  Fig.  1.  The  constitu¬ 
tive  relations  for  the  one-dimensional,  elastic,  anisotropic, 
micropolar  system  considered  in  this  work  are  given  by 

{T}„  =T  =  Ae+Cy,  {M}n  =  M  =  Cs+By.  (2.3) 

The  constitutive  parameters  A,  B,  and  C  describe  the  elastic  behav¬ 
ior  of  the  model  and  control  the  coupling,  or  interaction,  between 
the  two  kinematic  modes:  extension  and  rotation,  see  the  terms 
Ci  — C4  in  Shahsavari  and  Ostoja-Starzewski  (2005).  (The  tilde  over 
the  model  parameters  indicates  that  these  are  dimensional  quanti¬ 
ties;  non-dimensional  parameters  will  be  introduced  shortly.) 
Applying  Eqs.  (2.2)  and  (2.3)  to  Eq.  (2.1)  leads  to  the  following  sys¬ 
tem  of  PDEs: 


~S  2u  ^d2m  82u 
6x2  Sx2  F  0t2 


0^ 

0t2 


(2.4) 


where  for  a  one-dimensional  system  the  microinertia  term  jlk  from 
Eq.  (2.1 )  becomes  the  scalar  j.  Eq.  (2.4)  are  second-order,  linear  PDEs 
that  require  two  initial  conditions  and  two  boundary  conditions  for 
each  kinematic  quantity,  u  and  cp.  Eq.  (2.4)  may  be  written  in  the 
following  non-dimensional  form: 


AV2u  +  CV2cp  =  ii,  CV2u  +  BV2(p  =  Dtp, 
where 


B  = 


B  x2 

-ps*' 


(2.5) 


(2.6) 


The  terms  A,  B,  C,  and  D  are  non-dimensional;  x  is  unit  time  and  <5  is 
unit  length.  In  addition,  the  derivatives  indicated  by  V  and  the  dots 
are  now  derivatives  with  respect  to  position  normalized  by  S  and 
time  normalized  by  x,  respectively;  x  and  t  will  now  be  taken  as 
non-dimensional  quantities  normalized  by  S  and  x. 

One  form  of  the  strain  energy  density  that  yields  the  desired 
system  of  governing  equations  following  Hamilton's  principle  is  gi¬ 
ven  by 

W  =  WVu)2  +  2  CVuVcp  +  B(Vcp)2],  (2.7) 


2.1.  Micropolar  continuum  model 


such  that  the  (force)  stress  and  the  couple  stress  are  equivalent  to 
Eq.  (2.3)  and  are  given  by 


The  balance  of  linear  momentum  and  the  balance  of  angular 
momentum  equations  for  a  micropolar  body  are  given  by  Eringen 
(1999)  as  follows: 

02  .  02 

{T}w,k  —  P^p{U}l!  \M}kl,k+ elkm{T}km  —  PJlk^Wjk’  (21) 

where  the  second-order  tensors  T  and  M  and  the  vectors  u  and  tp  are 
defined  overx  e  Y\  t  e  T+.  The  micropolar  body  occupies  the  region 
W  in  Euclidean  space;  time  T+  is  defined  over  the  range  [0,  oo).  For  a 
micropolar  body,  there  are  two  stress  tensors:  the  (force)  stress  ten¬ 
sor  T  and  the  couple  stress  tensor  M.  In  addition,  there  are  two  vectors 
describing  the  deformations  of  a  micropolar  body:  the  displacement 
vector  u  and  the  microdisplacement  vector  <p.  The  mass  density  p  and 
the  microinertia  jlk  also  appear  in  Eq.  (2.1 ).  Standard  indicial  notation 
is  used,  i.e.,  {u}(  denotes  the  /th  component  of  the  vector  u,  the 
notation  k  denotes  the  partial  derivative  with  respect  to  position, 
and  £ikm  is  the  permutation  symbol.  Note  that  there  are  no  body  force 
or  body  couple  terms  included  in  Eq.  (2.1 ). 


dW 

~0fT 


M  = 


0W 

0y 


^(AVu  +  CV<p), 
^(CVu  +  BVp). 


(2.8) 


To  ensure  stability  of  the  thermodynamic  state  of  the  system,  it  is 
necessary  for  W  >  0  for  all  possible  applied  normalized  strains 
Vu  and  V<p,  see  Eringen  (1999).  The  following  inequalities  are  con¬ 
sequences  of  this  requirement: 

A  >  0,  B  >  0.  AB-C2>0.  (2.9) 


The  boundary  conditions  to  be  applied  to  Eq.  (2.5)  may  include  both 


1  For  the  one-dimensional  case,  there  is  no  distinction  in  the  governing  equations 
for  different  orientations  of  rotational  motion,  since  j,fe  =  0  for  /#/<  and  j//#0  for  only  a 
single  value  of  /.  There  will  always  be  two  governing  differential  equations  of  the  form 
of  Eq.  (2.4)  based  on  two  constitutive  relations  of  the  form  of  Eq.  (2.3),  although  the 
subscripts  may  differ. 
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- 

X 

Fig.  1.  Owing  kinematic  quantities  u  and  cp  by  tracking  a  point  from  its  reference 
configuration  to  a  deformed  configuration. 

essential  and  natural  conditions.  For  example,  if  u  and  (p  represent 
the  normalized  deformation  and  rotation  at  a  boundary  and  if  T  and 
M  represent  the  normalized  (force)  stress  and  couple  stress  at  a 
boundary,  then  the  essential  boundary  conditions  consist  of 

u|x=o,i  =  u,  (p\x= o,,  =q>,  te  [0,  oo),  (2.10) 

while  the  natural  boundary  conditions  consist  of 

^Vii|x=0ii  +  CV<p|x=0,i=T,  CV u|x=0ji  +  BV <p|„=0,i  =  M ,  te[0,oo), 

(2.11) 

recalling  Eq.  (2.8).  The  non-dimensional  stress  terms  T  and  M  are 
obtained  by  normalizing  by  pS2r~2  and  p<53 t~2,  respectively.  Initial 
conditions  specify  the  initial  velocities  in  terms  of  both  kinematic 
quantities,  u  and  cp.  In  other  words 

u|t=0  =  <P\t=0  =  co,  X6[0.1],  (2.12) 

where  v  is  the  normalized  initial  axial  velocity  and  cb  is  the  normal¬ 
ized  initial  rotational  velocity.  In  addition,  it  is  assumed  that  at  time 
f  =  0  the  body  is  stress-free,  recall  Eq.  (2.8),  i.e. 

AVu|t=0  +  CV<p|t=0  =  0,  CVu|t=0  +  BV<p|t=0  =  0,  xe[0,l]. 

(2.13) 

It  is  apparent  that  C  controls  the  coupling  between  the  translation 
and  rotation;  it  will  be  shown  that  A  is  related  to  the  elastic  modu¬ 
lus  for  the  classical,  non-polar  case. 

2.2.  Boundary  and  initial  conditions  for  the  two-layer  problem 

Consider  the  bar  consisting  of  two  different  micropolar  layers, 
layer  1  and  layer  2,  rigidly  fixed  at  x  =  1  as  shown  in  Fig.  2.  The 
properties  for  layer  1  are  given  by  A i,  B ,,  C,,  and  Dp,  for  layer  2, 
these  properties  are  given  by  A2,  B2,  C2,  and  D2.  If  the  first  body 
is  of  length  /,  where  0  <  /  <  1 ,  and  the  second  body  is  of  length 
1  -  /,  it  follows  from  Eq.  (2.5)  that  the  governing  equations  for 
the  two  bodies  are 

A]V2Ui  +  CiV2^,  =  iii,  C,V2t(i  -EfiiV2^,  =  D,<p,,  (2.14) 

A2V2u2  +  C2X72tp2  =  u2,  C2V2u2  +  B2V2tp2  =  D2ip2,  (2.15) 

where  u,  and  q> ,  are  defined  over  the  range  x  e  [0,/]  and  u2  and  tp2 
are  defined  over  the  range  x  g  [l,  1],  Recalling  Eqs.  (2.12)  and  (2.13), 
the  initial  conditions  for  a  stationary  body  with  zero  stress  state  are 
given  by 


U!|t=0  =  u2|t=0  =  o,  tp\  |t=0  =  <p2\t=0  =  0,  (2.16) 

and 

^iV[q|t=0  +  Ci  V  qc,  |r=0  =  0,  C1Vu1|t=0-|-B1V(p1|t=0  =  0,  (2.17) 

A2  Vu2|t=0  +  C2V<p2|t=0  =  0,  C2Vu2|t=0  +  B2V<p2|t=0  =  0.  (218) 

Instead  of  using  Eqs.  (2.17)  and  (2.18),  but  consistent  with  those 
equations,  we  will  assume  that  the  bodies  are  initially  undeformed, 
i.e. 

Ui|t=o  =  u2|t=o  =  0,  tPilt-0  =  <P2\t-0  =  0-  (2.19) 

The  essential  boundary  conditions  at  x  =  1  are  given  by 
u2|x=1  =  0,  <p2|x=1  =  0,  (2.20) 

while  the  natural  boundary  conditions  at  x  =  0  are  given  by 
^iVu1|x=0  +  CiV(p1|x=0  =  T,  C,Vu,|x=0 +  B1V<p1|x=0  =  M, 

(2.21) 

recalling  Eqs.  (2.10)  and  (2.11).  Finally,  the  matching  boundary  con¬ 


ditions  at  the  interface  x  =  l  are  given  by 

ui|x=,  =  u2|x=„  (2.22) 

<7>i  lx=i  =  (P2\x=i,  (2-23) 

AiViq^  +  CiV^U  =A2Vu2|x=(  +  C2Vr/>2|x=J,  (2.24) 

C,  Vu,|x=l  +  B,  Vtp,  |x=l  =  C2Vu2|x=,  +  B2  V  qo2  |X=J,  (2.25) 


for  t  g  [O.oo).  The  complete  mixed  boundary-initial  value  problem 
is  given  by  Eqs.  (2.14)-(2.16),  (2.20)-(2.25). 

2.3.  Boundary  and  initial  conditions  for  the  impact  problem 

Consider  now  the  impact  between  two  micropolar  bodies,  as 
shown  in  Fig.  3.  This  problem  is  similar  to  that  presented  in  Section 
2.2  in  that  there  are  two  different  bodies  with  two  sets  of  material 
properties  given  by  the  constants  Au  B1t  Cj,  and  for  the  initially 
moving  body  (the  flyer)  and  A2,  B2,  C2,  and  D2  for  the  initially  sta¬ 
tionary  body  (the  target).  The  governing  equations  are  given  by 
Eqs.  (2.14)  and  (2.15). 

The  first  body  is  traveling  at  an  initial  normalized  velocity  v 
and,  at  f  =  0,  it  impacts  the  target.  Therefore,  the  initial  conditions 
must  specify  both  the  initial  velocities  as  in  YuFeng  and  DeChao 
(1998)  and  Goldsmith  (1999): 

t<i  |t=0  =  vx  =  v,  (p,  |t=0  =  ffl,  =  0,  (2.26) 

“2|t=o  =  f>2  =  0.  <p2|t=0  =  ®2  =  0,  (2.27) 

as  well  as  the  initial  undeformed  state,  which  is  described  by  Eq. 
(2.19).  As  shown  in  Fig.  3,  the  free  ends  of  both  bodies  (x  =  0  and 
x  =  1 )  are  stress-free  leading  to  the  boundary  conditions 

A,Vu1|x=0  +  C1V(p1|x=0  =  0.  C1Vti1|x=0  +  B1V<p1|x=0  =  0.  (2.28) 
712Vu2|x=1  +  C2V<p2|x=1  =  0,  C2Vu2|x=1 +B2V(p2|x=1  =  0.  (2.29) 

The  matching  conditions  at  x  =  /  are  identical  to  those  given  by  Eqs. 
(2.22)-(2.24)  and  (2.25).  The  mixed  boundary-initial  value  problem 


T,  M 


f 

e 

< 8 

- 

(  /  j 

L 

_ ! _ J 

— nr 

Fig.  2.  A  body  consisting  of  two  different  micropolar  layers  is  subjected  to  a  known 
normalized  (force)  stress  and  couple  stress  at  x  =  0  and  is  rigidly  fixed  at  x  -  1 . 


t  <  0 
v 


h—i-H  h — i-i— H 


t  =  o 
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Fig.  3.  Two  micropolar  bodies  are  shown  before  impact  (t  <  0)  and  at  the  moment 
of  impact  (t  =  0).  The  first  body  of  length  l  (the  flyer)  is  traveling  with  a  normalized 
velocity  v\  the  second  body  of  length  1  -  /  (the  target)  is  initially  at  rest. 
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is  thus  given  by  Eqs.  (2.14),  (2.15),  (2.19),  (and),  (2.23)-(2.27)  2.28  3.  Solving  the  boundary  value  problems  and  the  discrete 

(and  2.29).  problem 


2.4.  The  discrete  model 


Although  the  primary  goal  of  the  present  work  is  to  study  con¬ 
tinuous  problems,  a  brief  outline  of  an  equivalent  discrete  system 
may  be  helpful  to  provide  some  physical  interpretation  of  the  re¬ 
sults.  A  discrete  system  consisting  of  a  series  of  masses  and  con¬ 
necting  springs  is  shown  in  Fig.  4.  This  model  is  similar  to  that 
introduced  in  the  study  of  dynamic  behavior  of  granular  media 
by  Lisina  et  al.  (2001)  and  Pavlov  et  al.  (2006).  Each  mass  is  con¬ 
nected  to  each  neighbor  with  two  springs:  one  spring  with  spring 
constant  k0  connects  the  centers  of  each  mass  and  one  spring  with 
spring  constant  k ,  connects  the  corners  of  each  mass.  A  given  mass 
i  is  allowed  to  translate  horizontally  (u,)  and  rotate  in  the  plane  of 
the  page  about  its  center  (<p,).  In  future  work,  it  may  be  interesting 
to  consider  non-linear  springs,  e.g.,  non-monotone  springs  as  in 
Balk  et  al.  (2001 )  or  bi-stable  springs  as  in  Puglisi  and  Truskinovsky 
(2000)  or  Slepyan  et  al.  (2005). 

For  a  system  with  n  masses,  the  total  energy  of  the  system  is 
found  by  summing  the  energy  associated  with  the  deformation 
of  each  spring: 

O  =  ]T  ^k0(Ui+i  -  Uf)2  +  \k,  [ui+1  -  Ui  +  Ufa  -  <pM)}2,  (2.30) 

i=i  z  z 


where  only  linear  terms  are  retained  from  the  development  of  the 
equations  based  on  Fig.  4.  The  2 n  Lagrange's  equations  are  obtained 
by  taking  derivatives  of  U  with  respect  to  the  n  discrete  axial  dis¬ 
placements  and  the  n  discrete  rotations  corresponding  to  each  of 
the  n  masses. 

In  addition  to  obtaining  a  system  of  Lagrange’s  equations,  by 
applying  a  Taylor  series  expansion  of  the  displacement  terms  we 
may  also  obtain  the  continuum  version  of  the  discrete  model  in 
the  following  form: 


[l2  (k0  +  ki)] 


cPu 

0X2 


(— a/2/<  i) 


d2cp 

0x2 


=  m 


dhi 

8t2 


(-a/2/t,) 


0^ 

0X2 


+  (a2/2k,) 


d2cp 

0X2 


0> 

0t2  ’ 


(2.31) 

(2.32) 


where  m  and  /  refer  to  the  mass  and  the  mass  moment  of  inertia  of 
each  mass  in  the  system.  Comparing  Eqs.  (2.31)  and  (2.32)  with 
(2.4),  it  follows  that  the  constitutive  parameters  from  the  contin¬ 
uum  case  may  be  related  to  the  discrete  parameters  as  follows: 

A  =  l2(k0  +  k,)/V,  B  =  a2l2kj/V,  C  =  ~al2ki/V.  p  =  m/V.  j  =  I/m , 

(2.33) 


where  V  is  the  volume  of  the  continuum  system  obtained  from  the 
discrete  system.  Since  k0,ki,a,  I,  V  >  0,  the  inequalities  from  Eq. 
(2.9)  are  identically  satisfied  by  Eq.  (2.33). 


fci 


Fig.  4.  Three  discrete  masses,  i  —  1 ,  i,  and  i  +  1,  each  allowed  to  translate  axially  and 
rotate  about  the  mass  center,  are  shown.  One  spring  (k0)  connects  the  centers  of 
adjacent  masses  together  and  another  spring  (ft, )  connects  opposing  corners  of 
adjacent  masses. 


In  this  section,  three  different  solution  techniques  are  described 
for  the  different  problems  presented  in  Section  2.  In  Section  3.1,  the 
D’Alembert  method  is  used  to  solve  the  bar  problem  described  in 
Section  2.2.  This  approach  is  particularly  useful  in  obtaining  ratios 
of  incident,  reflected,  and  transmitted  stress  waves  (reflection  and 
transmission  coefficients)  across  interfaces  separating  materials 
with  different  impedance  values.  A  Laplace  transform  solution  is 
described  for  the  impact  problem  in  Section  3.2  and  a  brief  descrip¬ 
tion  of  the  solution  method  for  the  discrete  problem  of  Section  2.4 
is  included  in  Section  3.3. 


3.1.  Solving  the  two-layer  problem  via  the  method  of  D'Alembert 


One  approach  that  may  be  used  to  solve  the  boundary  value 
problem  from  Section  2.2  is  to  use  a  D’Alembert  formulation. 
Recalling  Fig.  2,  solutions  for  Eqs.  (2.14)  and  (2.15)  for  the  case 
where  both  layers  in  a  bar  begin  at  rest  are  assumed  to  be  of  the 
form 

ui(x,  t)  =  px[Fi(t/Af  +  x)  +  F2(t/Af  ~x)\  +  p2[F3{t/Af  +x) 

+  F4(t/4-x)],  (3.1) 

<Pi(x,t)  =  F\[t/Xf  +x)  +  F2(t/A/  -x)  +  F3(t/Af  +x) 

+  F4(t/Aj  -  x),  (3.2) 

u2(x, t)  =  p3[F5(t/A2  +x)  +F6{t/A2  -x)]  +  p4[F7(t/A$  +x) 

+  F8(t/4-x)  ],  (3.3) 


cp2(x,t)  =  F5(t/A2  +x)  +F6(t/A2  —  x)  +  F7(t/ A2  -F x) 

+  Fs(t/4-x).  (3.4) 

The  following  constants  are  used  in  writing  Eqs.  (3.1 )— (3.4): 


Mi 


(2,)  -  d4  _(2+)  -  d4  _(A2)  —  d8  _(A2)  -d. 


,  M2 


,M3 


d  7 


Ma 


(3.5) 


where 


A±  =  \jdi+dA±  v/d2  +  4d2d3-2d1d4  +  d2/%/2,  (3.6) 

A2  =  '\J 1  ds  +  d8  ±  \J d2  +  4d6d7  —  2d$dg  +  djj/ \/2,  (3-7) 

and  the  constants  di,d2, _ ds  are  functions  of  the  model  parame¬ 

ters  as  follows: 


d,  = 


d4  — 


d7  = 


B , 


AiBj  -  C2 
A,Di 

AiBi  -  C2 
-C2 


,  d2  = 
,  d5  = 

,  ds  = 


-CiDi 

d  -c' 

(3.8) 

AiBi-C]' 

AiBi-C 2’ 

b2 

-C2D2 

(3.9) 

a2b2  -  C2  ’ 

D  A2B2  -  c2 ' 

a2d2 

(3.10) 

a2b2  -  c2 ' 

Based  on  the  forms  of  Eqs.  (3.6)  and  (3.7),  it  follows  that  A+>A,  and 
A2  >  A/ .  All  eigenvalues  are  real  as  a  consequence  of  the  inequalities 
given  by  Eq.  (2.9).  The  normalized  velocities  of  the  fast  waves  are 
1  /A/  and  1  /A2  and  the  normalized  velocities  of  the  slow  waves 
are  1  /Af  and  1  /Af.  The  ratios  of  the  fast  wave  speeds  to  the  slow 
wave  speeds  are  thus  Af  /A/  and  A2 /A/.  This  ratio  will  be  designated 
r  in  the  appendix,  see  Eq.  (A.l). 

For  a  non-polar  material,  there  is  infinite  rotational  stiffness  so 
that  there  can  be  no  rotation,  i.e.,  B  — >  oo.  In  this  situation, 
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1  /k~  — >  0  and  1  /X+  — >  VA.  To  write  the  wave  speed  in  dimensional 
form,  it  is  necessary  to  multiply  by  8/x,  which  is  unit  length  per 
unit  time.  Recalling  the  first  of  Eq.  (2.6),  we  see  that 


where  E  is  the  elastic  modulus,  p  is  the  mass  density,  and  \/Efp  is 
the  classical  wave  speed  for  a  one-dimensional  elastic  body.  This 
shows  that  in  the  non-polar  case,  A  corresponds  to  E. 

We  will  now  consider  a  simpler  case  with  identical  materials, 
i.e.,  the  bar  is  homogeneous  and  there  is  no  need  for  model  param¬ 
eter  subscripts.  It  is  then  only  necessary  to  solve  for  u(x,  t)  and 
tp(x,  t)  with  stress  applied  at  x  =  0  and  the  bar  rigidly  fixed  at 
x  =  1.  The  next  step  is  to  use  the  definition  of  (force)  stress  and 
couple  stress,  Eq.  (2.8),  with  the  D’Alembert  form  of  the  displace¬ 
ment  functions  and  the  boundary  conditions  Eqs.  (2.21)  and 
(2.20)  to  obtain 

(C  +  p,A)  VF,  (t/r )  -  (C  +  n,A) VF2(t/r)  +  (C  +  fi2A)VF3(t/X+) 
-(C  +  p2A)S7F4(t/X+)  =  f,  (3.12) 

(B  +  p,C)VFi (t/r)  -  (B  +  jutC)VF2(t/r)  +  (B+  fi2C)VF3(t/A+) 

-  (B  +  p2C)VF4(t/A+)  =  M,  (3.13) 

Pi F,  (t/A  + 1 )  +  /(, F2 {t/A~  - 1 )  +  fl2F3 ( t/X+  + 1 )  +  fl2F4(t/X+  - 1 )  =  0, 

(3.14) 

F,(t/r  +  1)  +F2(t/r  -  1)  +  F3(t/X+  +  1)  +  F4(t/X+  -  1)  =  0. 

(3.15) 

After  taking  the  Laplace  transform  of  Eqs.  3.12,  3.13  3.14,  and  3.15, 
we  may  solve  for  the  functions  Fi-F4  to  obtain 

oo 

F,  =  K\  ^{(2  +  4n  -  x  -  t/A~)H[-2  -  4n  +  x  +  t//T] 

n=0 

+  (-4  -  4n  +  x  +  t/r)H[- 4  -  4n  +  x  +  t/r]},  (3.16) 

OO 

F2  =  K\  J^{(2  +  4n  +  x  -  t/X~)H[-2  -  4n  -  x  +  t/r] 

n=o 

+  (-4n  -  x  +  t/r )H[-4n  -  x  +  t/r]},  (3.17) 

oo 

F3  =  K2  ^{(2  +  4n  -  x  -  t/2+ )H[- 2  -  4n  +  x  +  t/X+\ 

n= 0 

+  (-4  -  4n  +  x  +  t/A+)H[-4  -  4n  +  x  +  t/A+]},  (3.18) 

oo 


F4  =  k2  ]T{(2  +  4 n+x-  t/r )W [ — 2  -  4n  -  x  + t/X+] 


n=0 

+  (-4n  -  x  +  t/X+)FI[-4n  -  x  +  t/7+]}, 

(3.19) 

where 

A,/t2M  -B,T  +  C,(M  -  p2f) 

{Pi-p2){AiBi-C2i) 

(3.20) 

AiPiM  -  B]T  +  C^M  -  /(,T) 

(Hi-H2)(AiBi-C2i) 

(3.21) 

and  H[  ]  designates  the  Heaviside  unit  step  function.  The  D’Alembert 
form  of  the  solution  will  be  used  in  Section  4.5  to  obtain  ratios  of  inci¬ 
dent,  reflected,  and  transmitted  stress  through  an  interface  between 
materials  with  an  impedance  mismatch.  Although  the  summations  in 
Eqs.  (3.16)-(3.18)  and  (3.19)  are  infinite,  when  a  finite  time  t  is  con¬ 
sidered,  it  is  only  necessary  to  include  a  finite  number  of  terms  in 
the  summation  due  to  the  presence  of  the  Heaviside  functions. 


3.2.  Solving  the  impact  problem  via  the  Laplace  transform 

Laplace  transforms  will  be  used  to  solve  the  mixed  boundary- 
initial  value  problems  from  Sections  2.2  and  2.3,  although  the 
impact  example  from  Section  2.3  will  be  used  to  demonstrate  the 
method.  After  taking  the  transforms  of  Eqs.  (2.14)  and  (2.15)  and 
applying  the  initial  conditions  as  given  by  Eqs.  2.19,  2.26,  and 
2.27,  we  obtain  the  following  ordinary  differential  equations: 

AiX72u3  +C,V2^i  =s2ti1  -  y,C,V2Ui  +BiV2^>i  =  D1s2^>1,  (3.22) 
A2  V2u2  +  C2S72(p2  =  s2u2,C2'^2u2  +  B2V2(p2  =  D2s2cp2,  (3.23) 

where  the  bar  indicates  transformed  functions  and  s  is  the  variable 
in  the  Laplace  domain.  The  general  forms  of  the  solutions  to  Eqs. 
(3.22)  and  (3.23)  are 

Ui  =  +  c2/(1es;7x  +  c3^2e~s^*  +  c4pi2esiix  +  v/s2,  (3.24) 

(pi  =  Cie~sV  +  c2esri*  +  c3e~si(x  +  c4es>nx,  (3.25) 

u2  =  c5p3e-sl2x  +  c6/(3es,!2x  +  c7/j,4e^si2x  +  csfi4esl2x,  (3.26) 

(p2  =  c3  e~si2  *  +  c6  es22 x  +  c7  e~%x2 *  +  c8  x ,  (3.27) 

where  the  eigenvalues  and  constants  shown  are  defined  by  Eqs. 
(3.5)-(3.9)  and  (3.10).  The  eight  constants  ci,...,c8  are  found  by 
applying  the  general  solutions  given  by  Eqs.  (3.24)-(3.26)  and 
(3.27)  to  the  eight  boundary  conditions  given  by  Eqs.  (2.22)- 
(2.24),  (2.28),  (and)  (2.29).  One  consequence  of  applying  these 
boundary  conditions  is  that 

Ci  =  c2,  c3  =  c4,  c5  =  e2s27c6,  c7  =  e2si2c8.  (3.28) 

At  this  point  it  is  necessary  to  numerically  invert  the  solutions  from 
the  Laplace  domain  to  obtain  the  final  solution  in  the  time  domain. 
The  approach  used  in  this  work  is  the  Dubner-Abate-Crump  (DAC) 
algorithm  described  by  Crump  (1976);  the  effects  of  Gibbs  phenom¬ 
ena  in  the  solution  will  be  reduced  through  the  use  of  Lanczos’s  a- 
factors,  see  Laverty  and  Gazonas  (2006).  Results  from  the  Laplace 
transform  solution  are  used  in  Sections  4.2,  4.3  and  4.4. 

3.3.  Solving  the  discrete  problem 

Using  the  commercial  software  package  Mathematica,  an  expli¬ 
cit  formulation  is  employed  to  solve  the  discrete  problem.  First, 
displacements  u,  and  cpi  are  calculated  by  multiplying  the  longitu¬ 
dinal  velocity  Vi  and  rotational  velocity  co,  by  the  time  step  At  for 
each  mass  i.  Accelerations  are  found  by  determining  the  net  force 
and  moment  acting  on  each  mass  and  dividing  these  terms  by 
the  mass  m  and  mass  moment  of  inertia  /,  respectively.  The  forces 
and  moments  acting  on  each  mass  are  obtained  from  Eq.  (2.30). 
The  updated  velocities  V\  and  cuj  are  obtained  by  multiplying  each 
acceleration  by  At.  The  process  of  finding  displacements,  accelera¬ 
tions,  and  velocities  is  repeated  for  each  time  step  to  obtain  the 
dynamic  response  of  the  system.  Results  from  the  discrete  model 
are  used  in  Sections  4.1  and  4.2. 

4.  Results 

In  this  section,  five  different  examples  will  be  considered  that 
are  based  on  the  formulations  and  solutions  discussed  thus  far. 
The  first  example  in  Section  4.1  shows  a  discrete  model  and  it  is 
used  to  illustrate  the  presence  of  two  longitudinal  waves.  In  Sec¬ 
tion  4.2,  the  discrete  model  is  compared  with  an  equivalent  contin¬ 
uum  model.  An  example  of  impact  between  two  micropolar  bodies 
with  specially  chosen  model  parameters  is  presented  in  Section 

4.3.  In  Section  4.4,  some  additional  impact  examples  are  presented, 
including  one  that  demonstrates  the  effect  of  impedance  mismatch 
between  target  and  flyer.  Finally,  the  consequences  of  impedance 
mismatch  are  analyzed  further  in  Section  4.5.  In  all  plots  shown, 
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time  (the  horizontal  axis)  has  been  normalized  such  that  when 
normalized  time  equals  one,  the  slow  stress  wave  has  travelled 
the  length  of  the  specimen.  With  the  exception  of  Fig.  8,  the  verti¬ 
cal  axes  are  normalized  such  that  the  overall  length,  the  maximum 
velocity,  or  the  maximum  stress  are  set  equal  to  one,  depending  on 
the  results  being  plotted. 

4.1.  A  discrete  model 

This  example  is  the  discrete  equivalent  to  that  described  in  Sec¬ 
tion  2.2.  A  discrete  rod  made  up  of  80  discrete  masses  (n  =  80,  see 
Fig.  4  and  Eq.  (2.30))  with  one  end  rigidly  fixed  will  be  considered, 
i.e.,  Ugo  =  0  and  <p80  =  0.  At  the  free  end,  a  constant  compressive 
force  of  0.1  N  is  instantaneously  applied  at  time  t  =  0  s.  In  addition, 
the  following  model  parameters  have  been  chosen:  a  =  19.5  mm, 
b  =  3.5  mm,  /  =  20  mm,  k0  =  20  N/m,  lq  =  10  N/m,  m  =  0.01  kg, 
and  /  =  1CT6  kg  m2.  The  entire  structure  is  then  1.58  m  long.  (The 
numbers  were  chosen  to  represent  those  of  a  physical  system  that 
could  be  actually  constructed.)  The  positions  of  each  of  the  80 
masses  as  a  function  of  time  are  shown  in  Fig.  5.  It  is  possible  to 
observe  two  distinct  waves  emanating  from  the  origin  in  Fig.  5; 
the  faster  wave  reaches  the  fixed  mass  at  a  normalized  time  of 
approximately  0.55,  while  the  slower  wave  reaches  the  fixed  mass 
at  1.0  (the  time  scale  has  been  normalized  by  the  time  it  takes  for 
the  slow  wave  to  travel  the  bar  length).  A  similar  type  of  plot  is 
presented  in  the  work  of  Balk  et  al.  (2001). 

4.2.  Comparison  between  a  discrete  and  a  continuous  model 

In  this  section,  a  discrete  example  described  in  Section  2.4  is 
compared  with  a  continuous  example  from  Section  2.2.  The  prob¬ 
lem  to  be  considered  is  the  rigidly  supported  bar  shown  in  Fig.  2.  In 
the  present  example,  the  bar  is  homogeneous  and  there  is  no  need 
to  distinguish  parameters  with  subscripts.  From  Appendix,  the  fol¬ 
lowing  constants  for  the  continuous  case  are  used: 

A  =  10.8,  B  =  1,  C  =  -1.  D  =  0.02743.  (4.1) 

Using  Eqs.  (2.6)  and  (2.33),  the  following  equivalent  discrete  case 
constants  are  used:  a  =  1  m,  1  =  1  m,  k0  =  9.8  N/m,  /q  =  1  N/m, 
m  =  1  kg,  and  /  =  0.02743  kg  m2.  To  obtain  the  linear  version  of 
the  discrete  model  that  corresponds  to  the  linear  continuous  model, 
it  is  necessary  for  b  — >  0.  A  plot  showing  the  normalized  stress  at  the 
midpoint  of  each  rod  is  given  in  Fig.  6.  The  two  distinct  waves  are 
visible,  noting  that  the  faster  wave  is  twice  as  fast  as  the  slower 
wave  since  r  =  2,  see  Eq.  (A.l ).  With  the  exception  of  the  oscillatory 
behavior  of  the  underdamped  discrete  model,  the  two  solutions  cor- 


0  1  2  3  4  5  6  7 
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Fig.  5.  The  normalized  axial  positions  for  each  of  80  discrete  masses,  connected  as 
described  by  Fig.  4,  are  shown  as  functions  of  normalized  time  after  a  constant  load 
is  applied  to  the  end  i  =  1 .  The  end  i  =  80  is  fixed  for  both  axial  deformation  and 
rotation. 


Fig.  6.  Normalized  stress  at  the  midpoint  of  a  discrete  rod  and  a  continuous  rod  is 
shown  as  a  function  of  normalized  time.  Since  stress  at  the  midpoint  is  shown,  the 
slow  stress  wave  arrives  at  a  normalized  time  of  0.5. 

respond  relatively  well  through  the  first  three  reflections  of  the 
stress  waves,  e.g.,  after  a  normalized  time  of  two  for  the  fast  wave. 
After  this  point,  the  solutions  begin  to  diverge.  This  is  due  to  the  fact 
that  the  boundary  conditions  between  the  two  systems  are  funda¬ 
mentally  different.  For  the  discrete  model,  a  boundary  is  simply  a 
mass  at  an  end  that  is  connected  to  only  a  single  neighbor,  while 
an  interior  mass  is  connected  to  two  neighbors.  In  moving  to  the 
continuum,  the  interior  structure  of  Fig.  4  is  treated  as  a  representa¬ 
tive  volume  element  and  it  is  used  to  describe  the  behavior  of  the 
entire  body.  Hence,  in  Fig.  6,  each  time  a  wave  reflects  from  a 
boundary,  the  results  of  the  two  models  increasingly  differ. 

4.3.  An  impact  problem  with  tuned  micropolar  bodies 

The  case  to  be  considered  here  is  an  impact  between  two  iden¬ 
tical  micropolar  bodies  —  the  flyer  is  of  length  1  =  1/3  and  the  tar¬ 
get  is  of  length  1  -  /  =  2/3,  recalling  Fig.  3.  The  model  parameters 
have  been  chosen,  or  tuned,  to  ensure  that  the  first  occurrence  of  a 
net  tensile  stress  appears  at  the  interface  between  the  bodies  at 
x=  1/3.  (In  a  non-polar  material  system,  the  initial  net  tensile 
stress  will  appear  at  x  =  2/3,  which  is  the  middle  of  the  target.) 
The  following  constants  from  Appendix  are  used: 

Aj,A2  =  10.8,  B,,B2  =  1.  C,,C2=±1,  D,,D2  =  0.02743.  (4.2) 

Fig.  7  is  a  shock  wave  position-versus-time,  or  x-t,  diagram  that 
shows  the  propagation  of  the  stress  waves  (the  faster  wave  as  a 


Fig.  7.  The  x-t  diagram  for  the  impact  example  with  tuned  parameters  such  that 
the  initial  net  tensile  stress  occurs  at  point  (e),  the  interface  between  flyer  and 
target. 
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solid  line  and  the  slower  wave  as  a  dashed  line)  through  both  the 
flyer  and  target  immediately  after  impact.  The  points  in  the  figure 
are  labeled  to  correspond  with  the  points  to  be  shown  in  Fig.  8. 
Point  (e)  denotes  the  time  when  a  tensile  stress  state  initiates  at 
the  interface  between  the  flyer  and  the  target.  Only  by  assuming 
that  the  flyer  and  target  remain  attached  can  the  stress  waves  con¬ 
tinue  to  travel  through  the  interface  as  they  are  shown  to  do  in 
Fig.  7. 

in  a  plate  impact  experiment,  e.g.,  Antoun  et  al.  (2003),  the 
velocity  of  the  free  end  of  a  target  plate  (at  x  =  1  for  the  example 
presented  here)  is  measured  as  a  function  of  time.  The  results  of 
such  an  experiment  are  used  to  characterize  the  spall  strength  of 
the  target  material.  Fig.  8  shows  such  an  example  for  a  micropolar 
target  and  flyer,  characterized  by  Eq.  (4.2).  If  the  target  and  flyer 
were  to  separate  when  a  tensile  stress  first  appears  at  the  interface 
(at  point  (e)  in  Fig.  7),  Fig.  8  would  only  be  valid  up  to  the  velocity 
jump  at  point  (c).  The  results  of  Fig.  8  suggest  that  a  micropolar 
body  of  the  type  described  in  Section  2.1  should  posses  two  dis¬ 
tinct  jumps  in  velocity,  at  points  (a)  and  (b),  that  are  associated 


Fig.  8.  The  normalized  velocity  at  the  point  x  =  1  as  a  function  of  normalized  time, 
assuming  that  the  two  bodies  remain  permanently  attached  from  the  moment  of 
impact.  The  points  (a)-(e)  are  the  same  as  those  shown  in  Fig.  7. 


with  the  distinct  waves  discussed  in  Section  4.1.  The  maximum 
velocity  does  not  reach  the  initial  velocity  since  the  impact  also 
produces  rotational  motion,  in  other  words,  the  initial  purely  axial 
motion  generates  both  axial  deformation  as  shown  in  Fig.  8  as  well 
as  rotational  motion.  If  there  were  no  rotational  motion,  as  in  the 
case  of  a  non-polar  body,  then  u.2/v  =  \  at  its  maximum. 

In  order  to  see  the  stress  wave  propagation  for  this  example  in 
greater  detail,  Fig.  9  shows  four  snapshots  in  time  of  the  stress 
state  in  the  flyer  and  the  target.  The  arrival  of  the  first  stress  wave 
at  x  =  1  (at  point  (a)  in  Figs.  7  and  8)  is  shown  in  the  second  of 
these  figures.  The  initiation  of  a  tensile  stress  state  at  the  interface 
is  seen  to  occur  at  the  same  time  as  the  arrival  of  an  initial  com¬ 
pressive  wave  (the  larger  wave)  and  a  reflected  tensile  wave  (the 
smaller  wave)  at  the  free  end,  x  =  1  (at  Points  (b)  and  (e)  in  Figs. 
7  and  8).  The  fact  that  the  stress  waves  begin  as  compressive  waves 
and  are  reflected  from  the  free  surfaces  x  =  0  and  x  =  1  as  tensile 
waves  is  apparent  from  Fig.  9  as  well. 

Although  analytical  solutions  involving  a  finite  number  of  unit 
step  functions  have  been  used  thus  far,  it  will  be  necessary  to  re¬ 
sort  to  numerical  inversions  of  the  Laplace  transforms  in  Section 
4.4.  Therefore,  one  final  example  is  shown  in  Fig.  10  comparing 
the  numerical  solution  to  the  exact  solution  that  was  presented 
in  Fig.  8.  The  numerical  solution  is  found  using  the  DAC  algorithm 
with  Lanczos’  er-factors  with  512  terms  and  a  tolerance  equal  to 
1CT4,  see  Laverty  and  Gazonas  (2006).  As  can  be  seen,  the  effect 
of  the  Gibbs  phenomenon  is  minimal  and  so  we  can  be  confident 
in  the  results  presented  in  the  following  section  that  make  use  of 
this  numerical  tool. 

4.4.  Impact  with  different  materials 

Since  the  plate  impact  test  mentioned  in  the  previous  section  is 
a  common  approach  in  the  study  of  impact,  results  of  the  type  pre¬ 
sented  in  Fig.  8  will  be  considered  in  the  examples  to  follow.  The 
normalized  velocities  at  the  free  end  (x  =  1 )  of  a  target  impacted 
by  a  flyer  for  three  different  cases  are  shown  in  Fig.  11.  The  three 
cases  include  the  following:  a  non-polar  flyer  and  non-polar  target 
(np-np),  a  micropolar  flyer  and  micropolar  target  (mp-mp),  and  a 
non-polar  flyer  and  micropolar  target  (np-mp).  The  micropolar 
material  is  described  by  Eq.  (4.2),  while  the  non-polar  material  is 
modeled  with  the  following  parameters: 
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Fig.  9.  The  normalized  stress  state  within  the  system  at  four  different  times,  beginning  with  the  top  left  figure  and  ending  with  the  bottom  right  figure,  is  shown  as  a  function 
of  normalized  position. 
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Fig.  10.  To  illustrate  the  effect  of  Gibbs  phenomenon,  the  exact  solution  and  the 
numerical  solution  from  the  DAC  algorithm  are  plotted  for  the  case  shown  in  Fig.  8. 
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Fig.  11.  Normalized  velocity  at  the  point  x  =  1  as  a  function  of  time  as  was  shown 
in  Fig.  8.  In  this  figure,  three  examples  are  shown:  non-polar  flyer— non-polar  target, 
micropolar  flyer— micropolar  target,  and  non-polar  flyer— micropolar  target. 

Ax  =  10.8,  B,  =  10s,  C,  =  ±1(T6,  D,  =  0.02743  x  10s.  (4.3) 

Compared  with  the  micropolar  case,  the  coupling  constant  Ci  is  re¬ 
duced  by  six  orders  of  magnitude  to  reduce  the  coupling  between 
axial  deformation  and  rotation.  In  addition,  the  inertia  and  stiffness 
associated  with  rotation  are  increased  by  six  order  of  magnitude.  In 
terms  of  its  ability  to  rotate,  the  body  becomes  extremely  stiff  (due 
to  the  increase  in  B , )  and  insensitive  to  applied  rotation  deforma¬ 
tions  (due  to  the  large  microinertia  that  increases  D i). 

According  to  Fig.  1 1 ,  the  initial  response  to  impact  for  both  micro- 
polar  targets  is  similar,  exhibiting  two  distinct  velocity  jumps  and 
identical  wave  speeds.  This  is  due  to  the  fact  that  the  target  materials 
are  identical.  The  mismatch  in  material  properties  between  flyer  and 
target  for  the  “np-mp"  case  leads  to  a  more  complex  response,  as 
waves  are  partially  reflected  and  partially  transmitted  at  the  inter¬ 


face  between  dissimilar  materials.  A  detailed  analysis  of  the  reflec¬ 
tion  and  transmission  of  waves  across  an  interface  with  a 
mismatch  in  impedance  is  presented  in  Section  4.5. 

4.5.  Transmitted  and  reflected  stress  at  the  interface  of  a  Two-Layer 
micropolar  body 

In  addition  to  simply  using  the  D’Alembert  approach  to  solve 
the  given  boundary  value  problem,  this  approach  also  gives  insight 
into  the  relationships  between  incident,  transmitted,  and  reflected 
stress  waves  resulting  from  an  impedance  mismatch  at  the  inter¬ 
face  of  two  different  materials  described  by  different  constitutive 
parameters.  The  interface  is  perpendicular  to  the  direction  of  travel 
of  the  plane  stress  waves.  Eringen  (1999)  considered  harmonic 
plane  waves  traveling  through  micropolar  bodies  that  reflect  off 
of  a  free  surface,  while  Ghosh  et  al.  (2001)  examined  such  waves 
transmitted  across  a  boundary  between  two  micropolar  bodies. 

For  example,  consider  the  case  of  two  semi-infinite  one-dimen¬ 
sional  micropolar  bars  with  an  interface  at  x  =  0,  see  Fig.  12(a)  or 
(b).  If  the  eight  functions  Fi-F8  are  known,  the  behavior  of  each 
layer  is  completely  described  using  Eqs.  (3.1)-(3.3)  and  (3.4).  If 
we  only  wish  to  solve  for  the  reflection  and  transmission  coeffi¬ 
cients,  it  is  not  necessary  to  obtain  all  eight  functions  Fj-F8.  For 
example,  since  both  bodies  are  semi-infinite,  there  is  no  need  to 
consider  boundaries  at  ±oo.  Only  the  matching  boundary  condi¬ 
tions  at  x  =  0  given  by  Eqs.  (2.22)-(2.24)  and  (2.25)  are  used. 

Based  on  the  discussion  in  Section  4.1,  there  is  a  fast  wave  and  a 
slow  wave  and  each  is  capable  of  traveling  towards  the  left  and  the 
right.  Therefore,  for  each  layer,  there  are  four  waves  and  four  F-terms, 
Fi-F4  for  Layer  1  andF5-F8  forlayer  2.  For  example,  if  the  fast  wave  in 
layer  1  traveling  to  the  right  is  the  incident  wave,  then  we  are  also  only 
concerned  with  the  reflected  fast  and  slow  waves  and  the  transmitted 
fast  and  slow  waves.  In  this  particular  case,  f2  corresponds  to  the  inci¬ 
dent  fast  wave,  Fi  corresponds  to  the  reflected  fast  wave,  F3  corre¬ 
sponds  to  the  reflected  slow  wave,  F6  corresponds  to  the 
transmitted  fast  wave,  and  Fs  corresponds  to  the  transmitted  slow 
wave,  see  Fig.  12(a).  The  terms  containing  F4,  F5,  and  F7  may  be  ne¬ 
glected  and  we  are  left  with  four  equations  for  five  unknowns.  Similar 
relationships  hold  for  a  slow  incident  stress  wave,  as  shown  see 
Fig.  1 2(b),  and  for  incident  waves  traveling  towards  the  left  in  layer  2. 

For  the  case  shown  in  Fig.  12(a)  and  recalling  Eqs.  (2.22)-(2.24) 
and  (2.25),  the  four  equations  to  be  solved  to  find  F ,,  F3,  F6,  and  Fs 
as  functions  of  F2  are 


^[VF,+VF2]+^V3  =  ^VF6+^VF8,  (4.4) 

Xj  X,  x2  x2 

^[VF,+VF2]+1^V3=^VF6  +  irVF8,  (4.5) 

Xj  X]  x2  x2 

Ajj^IVF,  -  VF2]  +  /t2VF3}  +  Cj[VF,  -  VF2  +  VF3] 

=  A2[-/t3 VF6  -  /t4VF8]  +  C2[-VF6  -  VF8],  (4.6) 

CiIaMVFi  -  VF2]  +  /f2VF3}  +  B1[VF1~VF2  +  VF3] 

=  Gl-^VFg  -  /(4VF8]  +  B2[-VF6  -  VF8],  (4.7) 


Fast,  F2 - ► 

-Slow,  Fg 


-Fast,  Fs 


x  =  0 
(a) 


Slow,  Fa - >■ 

-Slow,  Fs 


-Fast,  Fs 


x  — >  00  x  - 
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Fig.  12.  Fig.  12(a)  shows  an  incident  fast  plane  wave  traveling  towards  the  right  in  layer  1  and  the  resulting  transmitted  and  reflected  plane  waves.  Fig.  12(b)  shows  the  same 
for  a  slow  incident  wave. 
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where  VFi  =  VF](t/7./),  VF2  =  VF2(t// l^),  VF3  =  VF3(t/2|), 
VF6  =  V/*6 (f/ A2 ),  and  VF8  =  VF8(t/72)  since  the  interface  is  taken 
to  be  at  x  =  0.  Eqs.  (4.4)  and  (4.5)  are  taken  from  Eqs.  (2.22)  and 

(2.23) ,  where  the  gradient  appears  after  taking  the  derivative  with 
respect  to  time  from  Eqs.  (2.22)  and  (2.23)  and  substituting  Eqs. 
(3.1)-(3.3)  and  (3.4).  Eqs.  (4.6)  and  (4.7)  come  directly  from  Eqs. 

(2.24)  and  (2.25). 

Eqs.  (4.4)-(4.6)  and  (4.7)  are  used  to  solve  for  VF i,  VF3,  VF6t 
and  VFS  in  terms  of  VF2,  i2,  Ah -7*4.  Aii2,  #1,2.  and  C3  2.  The  ratio 
of  the  fast  reflected  (force)  stress  wave  to  the  fast  incident  (force) 
stress  wave  is 


(a)  and  to  the  right  of  the  center  (b)  of  each  bar.  The  results  shown 
in  Fig.  13  are  obtained  by  solving  the  mixed  boundary-initial  value 
problem  via  the  Laplace  transform  approach  and  using  the  DAC 
algorithm.  By  observation,  the  results  of  both  approaches  are  the 
same,  thus  confirming  the  solution  techniques  employed.  That  is, 
the  reflected  (dashed)  stress  shown  in  Fig.  13(a)  has  a  magnitude 
of  0.116;  both  transmitted  (dashed)  waves  in  Fig.  13(b)  have  mag¬ 
nitudes  of  0.442. 

5.  Summary 


/VFA  Aj/tj  +  Cj  VF] 

vav 


(4.8) 


according  to  Eq.  (2.8).  The  ratios  for  the  slow  reflected  wave  to  the 
incident  wave,  the  fast  transmitted  wave  to  the  incident  wave,  and 
the  slow  transmitted  wave  to  the  incident  wave  are 


(VF3\ 

A,At2  +  Ci  , 

rVF6\ 

1*203  +C2 

rVF8\ 

1  712/(4  +  C2 

IvfJ 

Ai/q  +  Ci  ’  1 

VVFj 

'71,0,+ Ci’  1 

VVFj 

Ai/i,  +  Ci 

(4.9) 

respectively.  For  the  case  of  the  same  four  ratios,  but  now  with  re¬ 
spect  to  the  incident  slow  wave,  one  would  make  use  of  the  follow¬ 
ing  four  quantities; 


Ai/t, +Ci  VF3  - 

rVF6\ 

1  712A<3  +  C2 

rvFs) 

1  7\2A(4  +  C2 

VVfJ 

A,ai2  +  Ci  ’  VF4  ' 

'AiAG  +  Ci’1 

< vfJ 

71iA<2  +  Ci 

(4.10) 

since  F4  corresponds  to  the  incident  slow  wave. 

As  an  example,  consider  two  bars  with  a  normalized  (force) 
stress  applied  at  T  at  x  =  0  and  rigidly  fixed  at  x  —  1 ,  see  Fig.  2. 
One  bar  is  non-polar,  described  by  the  material  properties  of  Eq. 
(4.3),  and  the  other  bar  consists  of  two  layers.  The  bar  is  non-polar 
between  0  ^  x  ^  1/2  and  micropolar  between  1/2  ^  x  ^  1.  For  the 
two-layer  bar,  the  slow  stress  wave  is  the  incident  wave  transmit¬ 
ted  through  the  non-polar  body  that  reaches  the  interface  between 
the  layers.  By  considering  the  third  and  fourth  quantities  in  Eq. 

(4.10),  it  is  possible  to  find  material  parameters  that  will  reduce 
the  stress  transmitted  across  the  interface.  By  equating  the  stress 
corresponding  to  the  fast  and  slow  transmitted  waves,  it  is  possible 
to  obtain  the  following  parameters: 


A2  =  10.8, B2  =  102,C2  =  ±27.57. D2  =  0.02  743  x  102.  (4.11) 

By  examining  the  transmission  and  reflection  coefficients  obtained 
via  the  D’Alembert  method,  it  is  discovered  that  the  ratio  of  the 
transmitted  fast  and  slow  waves  to  the  incident  wave  is  0.442; 
the  ratio  of  the  reflected  slow  wave  to  the  incident  wave  is  0.116. 
These  results  are  shown  in  Fig.  13  for  a  homogeneous  bar  and  a 
two-layer  bar  at  the  positions  indicated  to  the  left  of  the  center 


A  one-dimensional  model  of  a  linear,  anisotropic,  micropolar 
body  subjected  to  transient  loading  associated  with  impact  has 
been  presented.  In  particular,  this  model  has  been  solved  for  the 
case  of  a  two-layer  bar  subjected  to  a  known  stress  at  one  end 
and  rigidly  fixed  at  the  other  end.  Either  one  or  both  layers  may 
be  considered  as  micropolar  bodies,  each  of  which  is  described 
by  four  model  parameters.  The  model  has  also  been  solved  for  an 
impact  problem  consisting  of  a  micropolar  projectile  and  a  micro- 
polar  target.  These  results  may  be  used  to  model  high  strain  rate 
experiments,  such  as  Kolsky  bar  tests  and  plate  impact  tests,  or 
to  give  insight  into  how  such  materials  would  behave  as  armors 
subjected  to  blast  and/or  ballistic  impact.  A  discrete  model  has  also 
been  solved  and  the  model  parameters  used  in  the  micropolar  con¬ 
tinuum  model  have  been  related  to  the  properties  of  the  discrete 
model,  which  are  physically  more  understandable. 

By  analyzing  the  solutions  to  the  governing  PDEs,  we  have  been 
able  to  choose  model  parameters  to  control  the  reflection  and 
transmission  coefficients  resulting  from  an  impedance  mismatch 
at  the  interface  of  two  different  materials.  This  allows  us  to  control 
the  transmission  of  stress  that  results  from  blast  or  impact.  For  the 
case  of  impact,  we  have  also  been  able  to  control  where  the  tensile 
stress  wave  first  appears.  This  is  helpful  since  materials  often  fail  in 
tension  as  a  result  of  the  combination  of  reflected,  initially  com¬ 
pressive,  stress  waves  caused  by  impact.  These  preliminary  results 
suggest  the  possibility  of  designing  an  optimal  system  that  will 
best  withstand  the  high  strain  rate  loads  it  will  be  subjected  to 
in  service.  This  sort  of  optimization  has  been  done  for  the  case  of 
multi-layered  elastic  strips  by  Velo  and  Gazonas  (2003). 

The  analysis  presented  herein  has  been  limited  to  linear  behav¬ 
ior,  i.e.,  requiring  small  displacements  and  assuming  linear  rela¬ 
tions  between  stress  and  strain.  In  addition,  it  has  been 
necessary  to  assume  a  periodic  microstructure,  an  assumption  that 
may  be  valid  in  the  case  of  an  engineered  structure  but  less  realis¬ 
tic  for  the  case  of  a  pulverized  ceramic.  By  adding  heterogeneity  to 
the  discrete  model  presented  in  Section  2.4,  a  homogenization 
analysis  approach  may  be  used  to  examine  the  effects  of  a  less  reg¬ 
ular  structure,  although  a  fully  random  structure  necessitates  a  dif¬ 
ferent  approach,  see  Ostoja-Starzewski  and  Trqbicki  (1999), 
Ostoja-Starzewski  and  Trqbicki  (2003).  Even  with  these  approxi- 
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Fig.  13.  Normalized  stress  waves  at  x  =  0.4  and  x  =  0.6  are  shown  for  a  uniform,  non-polar  bar  and  a  two-layer  bar  (the  first  layer  is  non-polar  and  the  second  layer  is 
micropolar).  The  solid  lines  indicate  the  stress  state  at  x  =  0.4,  while  the  dashed  lines  indicate  the  stress  state  at  x  =  0.6.  The  arrows  indicate  the  directions  the  waves  are 
traveling. 


C.L.  Randow,  C.A.  Gazonas  /  International  Journal  of  Solids  and  Structures  46  (2009)  1218-1228 


1227 


mations,  the  current  work  serves  as  a  benchmark  for  moving  to 
higher  dimensions  and  more  complex  constitutive  relations,  where 
it  becomes  more  difficult  to  relate  micropolar  constitutive  param¬ 
eters  to  physical  systems.  In  this  way,  the  results  presented  here 
will  guide  the  development  of  more  advanced  models  better  capa¬ 
ble  of  describing  the  response  of  multi-scale  materials  to  impact 
and  blast  conditions. 

Acknowledgment 

The  first  author  acknowledges  that  this  research  was  supported 
in  part  by  an  appointment  to  the  U.S.  Army  Research  Laboratory 
Postdoctoral  Fellowship  Program  administered  by  the  Oak  Ridge 
Associated  Universities  through  a  contract  with  the  U.S.  Army 
Research  Laboratory. 

Appendix.  Choosing  constants  for  the  micropolar  model 

In  this  appendix,  the  parameters  given  in  Eqs.  (4.1)  and  (4.2) 
will  be  calculated.  The  goal  here  is  to  determine  model  parameters 
that  ensure  the  following: 

•  The  faster  wave  speed  will  be  twice  the  magnitude  of  the  slower 
wave  speed  (r  =  2,  see  the  following  justification  for  this 
requirement),  and 

•  the  stress  corresponding  to  the  slower  wave  speed  will  be  larger 
than  the  stress  corresponding  to  the  faster  wave  speed  for  the 
case  with  no  applied  couple  stress. 


Fast-Wave  Rotational  Particle  Speed  2+ 
Slow-Wave  Rotational  Particle  Speed  X~ 


(A.5) 


If  there  is  only  a  single  wave  present,  as  in  the  classical,  non-polar 
case,  then  A  — »  0  or  A  — >  oo.  By  applying  Eqs.  (3.5)  and  (A.2)  to  Eq. 
(A.3)  and  solving  for  D,  it  follows  that 


1  +TA 

r2  +  a 


(A.6) 


Based  on  the  inequalities  in  Eqs.  (2.9)  and  (A.1),  it  follows  that  C 
equals  zero  only  when  A  — >  0  or  A  — >  oo,  i.e.,  the  condition  of  a  single 
wave.  In  the  limit,  if  r  — >  oo,  C  equals  ±VAB  based  on  applying  Eq. 
(A.4)  to  (A.2);  this  is  a  violation  of  the  last  inequality  in  Eq.  (2.9). 
Therefore,  r  must  be  finite  and  greater  than  one.  Since 
A,B,T,  A  0,  C  must  be  real.  In  conclusion  then,  to  ensure  that 
C¥=  0,  it  is  necessary  for  0  <  A  <  oo  and  1  <  r  <  oo.  Now  we  may 
consider  actual  numbers  for  our  model.  In  addition  to  requiring 
r  =  2,  we  will  choose  A  =  0.1  so  that  TS|0W  >  Tfast.  Making  use  of 
Eqs.  (A.2)  and  (A.6),  it  follows  that 

C  =  ±  VAB  =  ±0.3043 VAB,  D  =  ®  =  (0.2963) | . 

(A.7) 

The  values  of  A  and  B  were  chosen  to  make  C  equal  ±1,  i.e., 
VAB  =  1 8 /-\/30,  such  that 

A=  182/30=  10.8.  B  =  1,  C  =  ±l,  0  =  20/729  =  0.02743. 

(A.8) 


In  a  traditional  plate  impact  test  involving  non-polar  materials,  the 
target  material  will  fail  under  a  tensile  load  generated  at  the  mid¬ 
dle  of  the  target  specimen,  assuming  that  the  target  is  twice  the 
thickness  of  the  flyer  and  both  materials  are  identical.  As  a  demon¬ 
stration  of  the  nature  of  the  micropolar  material  as  described 
herein,  we  have  decided  to  present  an  example  where  the  initial 
tensile  stress  appears  at  the  target-flyer  interface,  rather  than 
within  the  target.  Since  the  typical  target-length-to-flyer-length 
ratio  equals  two,  it  will  be  necessary  to  require  0  =  2.  There  are 
no  additional  limitations  imposed  on  the  solution  under  these 
requirements.  Since  the  two  bodies  are  identical,  there  is  no  need 
to  specify  different  properties  through  the  use  of  subscripts.  The 
ratio  of  the  fast  wave  speed  to  the  slow  wave  speed,  mentioned 
in  the  first  of  the  two  requirements  is  defined  as 

Fast  Wave  Speed  1/2“  ^ 

Slow  Wave  Speed  1  /2+  >  ’ 

since  2+  >  Jr,  see  Section  3.1.  Based  on  Eqs.  (3.6)— (3.9)  and  (3.10) 
and  the  requirement  on  r  in  Eq.  (A.l),  it  follows  that  the  model 
parameter  C  equals 

c  =  ±  sJab{ i  +  r4)  -  b2d-'t2  -  A2Dr 2.  (A.2) 

Based  on  the  first  of  Eq.  (2.8)  and  the  solutions  to  Eqs.  (2.14)  and 
(2.15)  for  identical  flyer  and  target  materials,  the  ratios  of  the  stres¬ 
ses  corresponding  to  the  fast  and  slow  waves  are  given  by 


Tfast 

+  c\ 

( ftC±  B\ 

ft 

Tslow 

\Ji2A  +  CJ 

Vft  C±  BJ 

ft’ 

Mfast 

Mslow 


=  1,  (A.3) 


for  the  case  with  no  applied  couple  stress.  Next,  the  ratios  of  the  par¬ 
ticle  speeds  of  the  two  axial  waves  with  no  applied  couple  stress  are 

Fast-Wave  Translational  Particle  Speed  jtx  f 2+\ 

Slow-Wave  Translational  Particle  Speed  ji2  \JT  J  ’ 

(A.4) 


These  model  parameters  are  rewritten  as  Eqs.  (4.1)  and  (4.2);  there 
are  no  units  associated  with  these  parameters  due  to  the  non- 
dimensional  formulation  of  Section  2.1. 
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NEWARK  DE  19716 


NO.  OF 
COPIES 

1 

1 

1 

2 


1 

1 

1 

1 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


UNIV  OF  DELAWARE 
CTR  FOR  COMP  ST  MATRLS 
J  GILLESPIE 
NEWARK  DE  19716 

COMPUTATIONAL  MECH 
CONSULTANTS 
J  A  ZUKAS 
PO  BOX  11314 

BALTIMORE  MD  21239-0314 


2  WASHINGTON  ST  UNIV 
INST  OF  SHOCK  PHYSICS 
Y  M  GUPTA 
J  ASAY 

PULLMAN  WA  99164-2814 

1  ARIZONA  STATE  UNIV 

MECHCL  &  ARSPC  ENGRG 
D  KRAJCINOVIC 
TEMPE  AZ  85287-6106 


LOUISIANA  STATE  UNIV 
R  LIPTON 

304  LOCKETT  HALL 
BATON  ROUGE  LA  70803-4918 

INST  OF  ADVANCED  TECH 
UNIV  OF  TX  AUSTIN 
S  BLESS 
H  FAIR 

3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

APPLIED  RSCH  ASSOCIATES 
D  E  GRADY 

4300  SAN  MATEO  BLVD  NE 
STE  A220 

ALBUQUERQUE  NM  87110 

INTERNATIONAL  RSRCH 
ASSOC  INC 
D  L  ORPHAL 
4450  BLACK  AVE 
PLEASANTON  CA  94566 

AKT  MISSION  RSRCH  CORP 
M  EL  RAHEB 
23052  ALCALDE  DR 
LAGUNA  HILLS  CA  92653 

UNIV  OF  ILLINOIS 

DEPT  OF  MECHL  SCI  &  ENGRG 

A  F  VAKAKIS 

1206  W  GREEN  ST  MC  244 

URBANA  CHAMPAIGN  IL  61801 

UNIV  OF  ILLINOIS 
ARSPC  ENGRG 
J  LAMBROS 

104  S  WRIGHT  ST  MC  236 
URBANA  CHAMPAIGN  IL  61801 


1  NORTHWESTERN  UNIV 

DEPT  OF  CIVIL  &  ENVIRON  ENGRG 
Z  BAZANT 

2145  SHERIDAN  RD  A135 
EVANSTON  IL  60208-3109 

1  UNIV  OF  DAYTON 
RSRCH  INST 

N  SBRAR 

300  COLLEGE  PARK 
MS  SPC  1911 
DAYTON  OH  45469 

2  TEXAS  A&M  UNIV 

DEPT  OF  GEOPHYSICS  MS  3115 
F  CHESTER 
T  GANGI 

COLLEGE  STATION  TX  778431 

1  UNIV  OF  SAN  DIEGO 

DEPT  OF  MATH  &  CMPTR  SCI 
AVELO 

5998  ALCALA  PARK 
SAN  DIEGO  CA  92110 

1  NATIONAL  INST  OF 

STANDARDS  &  TECHLGY 
BLDG  &  FIRE  RSRCH  LAB 
J  MAIN 

100  BUREAU  DR  MS  8611 
GAITHERSBURG  MD  20899-861 1 

1  MIT 

DEPT  ARNTCS  ASTRNTCS 
R  RADOVITZKY 
77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  MIT 

DEPT  MATLS  SCI  ENGRG 
E  THOMAS 

77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139 

2  MATERIALS  SCI  CORP 
A  CA1AZZO 

R  LAVERTY 

181  GIBRALTAR  RD 

HORSHAM  PA  19044 

2  DIR  USARL 
AMSRD  ARL  D 
C  CHABALOWSKI 
V  WEISS 

2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

ABERDEEN  PROVING  GROUND 

81  DIR  USARL 

AMSRD  ARL  WM 
SKARNA 
J  MCCAULEY 
P  PLOSTINS 
J  SMITH 
T  WRIGHT 
AMSRD  ARL  WM  B 
J  NEWILL 
M  ZOLTOSKI 
AMSRD  ARL  WM  BA 
D  LYON 

AMSRD  ARL  WM  BC 
P  WEINACHT 
AMSRD  ARL  WM  BD 
P  CONROY 
B  FORCH 

R  PESCE  RODRIGUEZ 
BRICE 

AMSRD  ARL  WM  BF 
W  OBERLE 
AMSRD  ARL  WM  M 
R  DOWDING 
S  MCKNIGHT 
AMSRD  ARL  WM  MA 
J  ANDZELM 
R  JENSEN 
A  RAWLETT 
M  VANLANDINGHAM 
E  WETZEL 

AMSRD  ARL  WM  MB 
M  BERMAN 
T  BOGETTI 


M  CHOWDHURY 
WDEROSSET 
W  DRYSDALE 
A  FRYDMAN 
D  HOPKINS 
L  KECSKES 
T  H  LI 

S  MATHAUDHU 
M  MINNICINO 
B  POWERS 
J  TZENG 

AMSRD  ARL  WM  MC 
R  BOSSOLI 
S  CORNELISON 
M  MAHER 
W  SPURGEON 
AMSRD  ARL  WM  MD 
J  ADAMS 
B  CHEESEMAN 
ECHIN 
KCHO 
B  DOOLEY 
C  FOUNTZOULAS 
G  GAZONAS 
J  LASALVIA 
P  PATEL 
C  RANDOW 
J  SANDS 
B  SCOTT 
CF  YEN 

AMSRD  ARL  WM  SG 
T  ROSENBERGER 
AMSRD  ARL  WM  T 
P  BAKER 
M  BURKINS 
AMSRD  ARL  WM  TA 
S  SCHOENFELD 
AMSRD  ARL  WM  TB 
N  ELDREDGE 
J  STARKENBERG 
AMSRD  ARL  WM  TC 
T  BJERKE 
T  FARRAND 
K  KIMSEY 
M  FERMEN  COKER 
D  SCHEFFLER 
S  SCHRAML 
S  SEGLETES 
AMSRD  ARL  WM  TD 
S  BILYK 
D  CASEM 
J  CLAYTON 
D  DANDEKAR 
N  GNIAZDOWSKI 


NO.  OF 

COPIES  ORGANIZATION 


M  GREENFIELD 
R  KRAFT 
B  LOVE 

M  RAFTENBERG 
E  RAPACKI 
M  SCHEIDLER 
T  WEERASOORIY A 
AMSRD  ARL  WM  TE 
J  POWELL 
B  RINGERS 
G  THOMSON 
AMSRD  ARL  WM  UV 
S  WILKERSON 
AMSRD  ARL  VT  RP 
J  BORNSTEIN 
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